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^ . Abstract 

. We study the dynamics of first-order cofinement-deconfinement phase transition through 
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sure tend to speed up the transition. Another interesting feature is the start of the 



nucleation of hadronic bubbles in an expanding quark gluon plasma in the context of 
heavy ion collisions for interacting quark and hadron gas and by incorporating the 
effects of curvature energy. We find that the interactions reduce the delay in the phase 
i transition whereas the curvature energy has a mixed behavior. In contrast to the case 

of early Universe phase transition, here lower values of surface tension increase the 
' supercooling and slow down the hadronization process. Higher values of bag pres- 
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\ hadronization process as soon as the QGP is created. 

1 Introduction 

The ongoing and planned experiments involving heavy ion collisions at RHIC, CERN's SPS 

^ ' and at LHC are aimed at studying the behavior of quantum chromodynamics (QCD) at 

high energies. It is believed that at these high energies [T] (100 GeV/nucleon at RHIC and 

3 TeV/nucleon at LHC in the cm. frame) a hot baryon free plasma of quarks and gluons 

(QGP) is expected to be created even as at lesser energies, nuclear collisions at Brookhaven's 

AGS and at SPS show large stopping of the nuclei and hint at the baryon rich matter. The 

QGP plasma will soon undergo a phase transition to hadronic matter. The phase transition 

from QGP to hadron resonance gas (HRG) has been extensively studied in the context of 

the early Universe where the time scale is of the order of 10"^ — 10"'' s [2-4]. Perturbative 

finite temperature QCD methods along with relativistic transport theory have been used to 

study the early stages of the nuclear collisions by computer simulations [5,6]. The initial 

quantitative picture that emerges is that of a QGP fireball of ~ 150/m^ created at the initial 
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temperature of ~ 300 — 350 MeV, which is about twice the expected critical temperature Tc. 
Studies of quark-hadron phase transition in the early Universe, in heavy ion collisions and 
in high density nuclear matter in the core of compact stars [Zj point to the importance of 
interactions in the two phases. It has been shown that unless hadronic interactions are taken 
into account, the hadronic phase again becomes thermodynamically the preferred phase at 
high temperature and that QCD interactions are important in the QGP phase [H] making 
the critical parameters depend on the interactions. The lattice results also point to the 
departure of the QGP equation of state (EOS) from the ideal gas EOS even at temperatures 
well above T^. However there still remain uncertainties in the equation of state in both the 
confined and the deconfined phases. It has also been pointed out that in situations in which 
nucleation takes place, the curvature energy 87r7r, a term in the free energy in addition 
to the surface energy term Anr'^a (where a and 7 are the surface and curvature energy 
densities respectively) also plays an important role and should be kept in the calculation of 
the nucleation rate [9-12]. It is in this context that the above experiments will go a long way 
in establishing the dynamics of the phase transition and the final picture that will emerge 
would be of great interest in astrophysics and cosmology. If the phase transition from quark 
gluon plasma to confined hadronic matter is a first order thermodynamic phase transition, 
the transition would proceed through the nucleation of hadronic bubbles. Conventionally, 
the QGP would supercool to temperature below only then the nucleation starts. The 
bubbles greater than the critical size would expand driving the QGP phase into smaller 
and smaller regions thereby converting more and more matter into the hadronic phase. The 
expansion of the bubble would be accompanied by the release of latent heat which will reheat 
the plasma shutting off further nucleation. If the phase transition has to be completed, the 
latent heat must be gotten rid of. Such studies have been made in the context of the early 
Universe where it has been shown that phase transition proceeds in the expanding Universe 
after the Universe supercools to a temperature below Tc, the expanding bubbles then reheat 
the Universe expelling the QGP phase till the entire Universe is converted to the hadronic 
phase. The time of completion and the degree of supercooling depends on the EOS and the 
parameters like the Bag pressure and the surface tension j^. 

Csernai and Kapusta [Hj have applied a recently computed nucleation rate to a first 
order phase transition in a set of rate equations to study the time evolution of QGP in 
heavy ion collisions. Based on Bjorken hydrodynamics and on current parameter values, 
they find the transition generates about 30 percent extra entropy and also a time delay 
of 11/m/c in completion of the transition. They also estimate that the system should 
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supercool about 20 percent below Tc before nucleation of hadronic bubbles is rapid enough 
to start reheating the system. Csernai and Kapusta for the purposes of illustration used 
noninteacting system of hadrons and QGP. They modelled the hadronic phase by a massless 
gas of pions and the plasma phase by a gas of gluons and massless quarks of two flavors. 
In the light of the importance of taking interactions both in the hadronic as well as in the 
plasma phase for purposes of investigating phase transitions at high energies, temperature 
and densities as discussed above, we calculate the rate of formation of the hadron bubbles 
for interacting QGP and hadron resonance gas. It is our endeavour here to examine closely 
the effects of interactions on the dynamics of phase transition as the hot plasma of quarks 
and gluons is created and the temperature drops to Tc where a phase mixture of QGP and 
HRG develops. For this purpose QGP is treated as a gas of massless u,d quarks, massive 
s quarks and massless gluons. QCD interactions are treated perturbatively to order g"^ 
and the long range confinement effects are parameterized by the bag pressure B. For the 
hadron gas we use the known spectrum of low lying baryons and meson resonances and the 
repulsive interactions between them by incorporating the hard core excluded volume effects. 
In the literature hadronic interactions have been taken in many different ways (see Ref.7 for 
example), hadronic interactions can be accounted for through an average mean field repulsive 
potential arising from the Reid potential and vr — vr effective interactions or alternatively, 
repulsive effects can also be included by modifying the chemical potential as done by Kapusta 
et al. in a simple method of implementing the mean field theory. One could also use the 
full fledged framework of relativistic mean field theory itself. It has however been shown by 
Goyal et al.[7] that in the baryon free case , all the interaction schemes mentioned above 
give roughly the same value of the hadronic pressure, there are however, marked deviations 
at very high T and/or at high fis regime. 

The effect of curvature energy term on the nucleation of hadronic bubbles has been stud- 
ied earlier [U]. It has been found that it leads to a minima in the free energy in addition to 
a maxima below the critical temperature and also a minima above Tc, giving rise to sponta- 
neous nucleation of equilibrium sized bubbles. This may contribute to the completion of the 
phase transition and has interesting consequences. We study here the effect of curvature on 
the dynamics of the phase transition for several values of the parameters and discuss their 
consequences. 

The plan of the paper is as follows. In section 2 we discuss the bubble nucleation along 
with the importance of the prefactor which plays a crucial role in guiding the behavior of the 
nucleation rate. It turns out that the exponential factor is of less relevance here unlike the 
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case of early Universe where exponent dominates the behavior of the nucleation rate. We 
also study here the effect of curvature on the nucleation rate and give the two rate equations 
to study the time evolution of the transition for different parameter values. In section 3 
we discuss our result and compare them with the case when interactions and the curvature 
energy are not present. We also look at the supercooling that takes place and evaluate the 
behavior of the system with different parameter values. Finally in section 4 we give our 
conclusions. 

2 Bubble nucleation 

When the QGP cools through the critical temperature Tc, the new hadron phase becomes 
preferred. Energetically the new phase remains unfavourable as there is surface free energy 
associated with the QGP and the hadron gas interface. Creation of bubbles of the new phase 
are thus unfavourable and all nucleated bubbles with radius less than the critical radius die 
out but those with radii grater than the critical radius expand untill they coalesce with each 
other. Thus supercooling occurs before the new hadron phase actually appears and takes 
over. Consequently reheating takes place due to the release of latent heat. The critical radius 
is obtained by extremizing the thermodynamic work expanded to create a bubble, i.e. 

W = -—r\Ph - Pg) + Anr'a - 87r(7, - 7,)r, (1) 

where Ph and Pg are the pressures in the hadronic and quark phase respectively, a is the 
surface tension and 7 = (7g — 7/t) is the curvature coefficient and has been estimated in 0. 
Following Mardor and Svetitsky the curvature energy in the MIT model for massless quarks 
is given by 

E,= ^ r dkk{l + exp{k - fiB)/T}-' (2) 
ovr Jo 

where g is the statistical weight and r the radius. We get 7 ~ 7^ = where we have 
assumed ■jh « 7g following jTTj and /i^ to be negligible compared to T. According to 
Madsen ^0] , the assumption of zero strange-quark mass in the curvature term should be 
relaxed but presently no calculation of the curvature term for massive quarks exists. The 
surface tension is mainly due to the finite strange quark mass as the massless particles do not 
contribute to the surface term in the free energy. The critical radius is obtained by putting 
^ = and we get an extremum below given as 

-c-, = -^(l + V^) (3) 
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where AP — Pf^ — and (3 — ^^7- Since (3 is positive definite, a real solution exists only 
if /3 < 1. This critical radius is actually a maximum. There is another extremum solution 
(a minima) given by the critical radius 

a 



(4) 



These bubbles correspond to equilibrium bubbles which stabilse to their critical size Vc- and 
do not grow. The free energies of the critical sized bubbles are given by 

^c+ = ^[2 + 2(l-/3p-3/?] (5) 

and 

^c- = ^,[^-'^i^-(^f^"-m (6) 

respectively for the two critical radii. (3 decides the effect of the curvature term and decreases 
the critical radius for the expanding bubbles. The existence of a critical radius corresponding 
to the minimum of the free energy is a consequence of the curvature term. In fact there is a 
extremum in the free energy even above the critical temperature for all starting temperatures 
of the plasma. This extremum is also a minimum and the corresponding critical radius and 
the free energy are given by 



^c+> = ^(-l + Vl + /5) (7) 

W^c+> = ^[2-2(l + /3p + 3/?] (8) 

Significantly, the restriction of /3 < 1 is not present for these bubbles. This feature of the 
curvature term ensures that the phase transition actually begins well above the critical tem- 
perature (as soon as the plasma is formed) by the nucleation of equilibrium sized bubbles of 
the hadron phase, and this continues till below the critical temperature when the expanding 
bubbles also join them and start the rapid nucleation and expansion of the hadron phase. 
As pointed out earher, below Tc if (3 becomes greater than one, nucleation of both types of 
critical bubbles dissapear. The phase transition can now take place by the expansion of all 
the already nucleated critical bubbles rc+, Tc- and rc+>. This is because the critical radius 
Tc and the critical free energy Wc necessary for creation of growing bubbles become zero. 
This essentially means that the phase transition now gets rolling and completes fast. Sooner 
13 becomes greater than 1, faster will this process dominate the phase tranisition. As we will 
see , for certain reasonable values of the parameters {B and a) (3 becomes greater than 1 
almost near Tc causing unrestrained expansion of all the bubbles (ones above Tc as well as 
whatever exists below it) and completion of the phase transition. 
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The bubble nucleation rate(number of bubbles formed per unit time per unit volume) at 
temperature T is given by 

/ = I.e-'^^/T (9) 

where Jq is the prefactor having dimensions of T^. The prefactor used traditionally in the 
early Universe studies is given by Jo = (VFc/27rT)3/2T^. Csernai and Kapusta |E] have 
recently calculated this in a course grain effective field theoretic approximation to QCD and 
give 

16 / a \^/^ o-rigTc 



where rjq = 14. 4T^ is the shear viscosity in the plasma phase, is a correlation length of 
order 0.7 fm in the plasma phase, and Aw is the difference in the enthalpy densities of the 
two phases. This nucleation rate is limited by the ability of dissipative process to carry 
latent heat away from the bubbles surface as indicated by the dependence on the viscosity. 
The above expression however gets modifed in the present case because of the presence of 
the equilibrium bubbles. We now have different J's depending on which kind of bubbles are 
being discussed, namely /+, /_ and /+> corresponding to critical radiuses rc+, rc_ and rc+> 
respectively. 

The pressure in the QGP phase can be calculated by using the thermodynamic potential 
given in Kapusta and Shuryak [7j . For the hadron phase we use the known masses of the 
low lying 33 baryons and 45 mesons whose masses and degeneracy factors are taken from the 
Particle Data Group Summary [13j. The hadronic pressure and number densities for non- 
interacting point baryons and mesons are given by the usual thermodynamic relations. One 
of the simplest ways to account for the short range repulsive forces as discussed above is by 
considering the finite volume of the baryons that modify the space available for occupation 
in a manner akin to that of a Van der Waal's equation. We take baryons and antibaryons to 
have the same size as protons, given hj Vp = rrip/AB where B is the bag pressure. A radius 
Vp lying between 0.6 fm and 0.8 fm is often used [7] . The hadronic pressure and baryon 
densities corrected for finite volume effects are now given by 

1 + Efe < 1 + ^5 Vp 



nB= . + . (11) 



where b, b and m stand for baryons, anti-baryons and mesons respectively. The pressure 
equilibrium between the two phases sets which is decided by the parameters in the theory. 
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essentially by the parameters B and a. With nucleation rate given by eq.(9), one would like 
to know the fraction of space h{t) which is in the hadronic phase at the given time t as 
measured in a local comoving frame of the expanding system. Guth and Weinberg J3] have 
proposed a formula for h{t) in the cosmological first order phase transition. Csernai and 
Kapusta [H] have given a more accurate kinetic equation than the one used by Guth and 
Weinberg and it allows for the transition to complete. If the system cools to at time t^, 
then at some later time t, the fraction of space which has been converted to hadronic gas is 
given by them as 

h{t) = f dt'I{T{t')){l - h{t')]V{t', t) (12) 



where V{t, t') is the volume of a bubble at time t which had been nucleated at the earlier time 
t'; taking bubble growth into account. Here we have neglected bubble collisions and fusion. 
The above expression however gets modifed in the present case because of the presence of the 
equilibrium bubbles. Now we have an additional contribution to the hadron fraction from 
the bubbles that are nucleated above the critical temperature plus those that are nucleated 
below the critical temperature and do not grow in size (unless (3 > 1). So, the new expression 
is now given by 

h{t) = t dt'I+{T{t')){l - h{t')}V{t\t) 

Jtc 

+ / dt'l4T{t')){l-h{t')}-rUnt')) 

Jtc O 

+ £ dt'iMnmi - hit')}fri^int')) (13) 

The dynamical equation which couples time evolution of temperature to h{t) is obtained by 
using the longitudinal scaling hydrodynamics of Bjorken. It is given by 

Tr-j 

where e(T) and w{T) are the energy and enthalpy densities respectively being given by 
e(T) = h{t)eh(t) + (1 — h(t))eg(t), Ch and being the hadron and quark energy densities 
at time t respectively, likewise for the enthalpy density. Following ^B] the growth velocity 
of bubbles is taken to be v(T) = vo{l — T/Tc)^^"^ where vq = 3c for T > 2/3Tc. This gives 
V < c/v^, the speed of sound for a massless gas. We have assumed that if supercooling 
takes T below 2/3Tc then v is given by c/ v^- This ensures that closer T is to slower the 
bubbles grow. The volume of an expanding bubble at time t is then given by 

V{t',t) = fir^m) + £ dt"v{T{t"))}' (15) 
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In calculations we use the parameters C,q = 0.7 fm and riq = 14AT^ as in a and B are 

varied to study the behavior of the phase transition as the QGP cools and hadron bubbles 

are formed. We solve the two coupled equations (13) and (14) for different values of the 

parameters B and a. The nucleation time is defined as 

-1 47r 3^ 
3 " 

neglecting bubble growth. In our case we have a contribution to r from /+, /_ below Tc and 
/+> above Tc. 



3 Results and Discussion 

We find that the inclusion of the interactions in both the plasma phase and the hadron gas 
phase lowers the critical temperature, so higher values of the bag pressure are required for 
having reasonable values of Tc. For example, the value i?^/^ = 300MeV gives a Tc ~ 188MeV 
whereas B^/^ = 235MeV gives a Tc ~ lAAMeV . We show below plots of free energy vs 
radius of a hadron bubble, nucleation time vs temperature, temperature vs time and hadron 
fraction vs time etc. In fig.l we plot the free energy of a hadron bubble against it's radius 
in the QGP plasma. For temperatures just above the critical temperature we see that 
the effect of the curvature term is to produce a minima in the free energy which is to be 
contrasted with no extremum without a curvature term. Fig. la pertains to T = 1.0 ITc 
and the dashed and the dotted curves are without the curvature term. It is clearly seen 
that equilibrium sized hadron bubbles appear even above the critical temperature (solid and 
long dashed curves). Fig. lb is plotted for a different set of parameters. The lower panel 
is plotted just below the critical temperature at T = 0.99Tc. Here we find that below the 
critical temperature the curvature effect produces a minima as well as a maxima in the free 
energy. The maxima bubbles grow whereas the minima bubbles reach an equilibrium size. 
The dashed and the dotted curves are without the curvature effect showing the presence of 
only the maximum sized bubbles. Figures Ic and Id plot the same curves but for a different 
set of the parameters. It is interesting to note that for 5^/^ = SOOMeV, a ~ IMeV/ fm? 
[a^f^ = 65MeV) and T = 0.99Tc the free energy has no extremums, only a maximum value 
at zero (long dashed), indicating that critical radius is zero and all existing bubbles can now 
expand to reduce the free energy. In such cases j3 becomes greater than 1 soon after Tc 
is crossed and there cannot be nucleation of any fresh bubbles. The expansion of all the 
existing bubbles ensures completion of the transition. In Fig. 2 we have shown the plot of 
T/Tc agains time. The temperature of the QGP plasma decreases with time till Tc, when the 
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phase transition conventionally begins and supercooling starts. The degree of supercooling 
depends on the parameters and can be a significant fraction of Tf.. The important feature to 
see here is that the prefactor Jq is the deciding factor for supercooling and the duration of 
the transition rather than the exponent. On the other hand in the early Universe studies, 
it is the exponent which decides the behavior of the transition. In fig. 2a the dashed curve 
(without curvature) and the solid curve (with curvature) are for parameters B^^^ = 235MeV 
and a = bQMeV / fw? {a^/^ ~ 12bMeV). As can be seen they overlap, indicating that in 
this case the curvature term does not make any difference. This value of the surface tension 
is rather large and in fact exceeds the upper limit set in [15^. However, for the parameters 
B^l'^ = 235MeV and a = 20MeV/fm^ (a^/^ ~ 92MeV) there is a delay in the completion 
of the phase transition when the 7 term is included (long dashed) compared to the no 7 
(dotted) case. The delay is of the order of 13 fm and is because of (3 becoming greater than 
one a little below Tc with consequent stoppage of all bubble nucleation. When the 7 term 
is not taken, there is no such stoppage and continous nucleation ensures faster completion. 
Figure 2b shows the same curves for the parameters B^^^ = 235MeV, a ~ 7MeV/fm^ and 
= 300MeV, a ~ 7MeV/ fm^. In the latter case, the curvature term completly inhibits 
bubble nucleation below Tc and only the equilibrium bubbles formed above Tc expand to 
complete the phase transition (long dashed). This corresponds to the case of no extremum 
in free energy (fig. 2b). The inclusion of 7 term in this case actually reduces the time to 
complete the phase transition. This can be contrasted with the case corresponding to the 
other set of parameters in fig. Id and the cases shown in fig.2a. Comparing with the no 
curvature case. It is seen that increasing B speeds up the phase transition whearas decreasing 
it slows it down if the 7 term is included. We also see that the decrease in the surface tension 
increases the supercooling of the plasma by significant amount (can be close to 50 percent). 
However, increasing B reduces supercooling when 7 is included. 

Next we look at the nucleation rate vs T/Tc. One can see from these graphs (figures 3a 
and 3b) that nucleation of hadronic bubbles can start even above the critical temperature 
when we include the curvature energy. The effect of decreasing a is to suppress nucleation. 
We may note that in the early Universe reverse is the case, and it happens because there the 
exponent drives the nucleation rate whearas in QGP the prefactor drives the nucleation. The 
point at which (3 becomes greater than one is clearly visible when curvature term is included. 
In figure 3a the long dashed curve terminates at the point where further nucleation stops. It 
is also clear that nucleation starts the moment the plasma is created well above the critical 
temperature if the curvature term is included. Figure 3b shows that for smaller surface 
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tension the j3 term becomes greater than one soon after crossing through Tc (sohd and the 
long dashed curves). The other two curves shown are without the curvature terms and so 
begin below the critical temperature. 

In figure 4 upper panel the critical radius of a bubble is plotted as a function of T /T^ 
for those bubbles which are nucleated at the minima in the free energy. These bubbles 
in general are nucleated above Tc(rc+>) as well as below Tc(rc_). In fig.4a we see that 
for the parameters B^^^ = 235MeV, a = 20MeV/ fm^ the bubbles keep nucleating from 
above the critical temperature to well below it. However for another set (B^^'^ = 300MeV, 
7MeV/fm^) these bubbles are no longer nucleated soon after crossing T^. In figure 4b 
we find this feature in both the parameter sets. In all such situations (where f3 becomes 
greater than 1) the completion of the phase transition is ensured by expansion of already 
nucleated bubbles. The lower panel plots the critical bubble radius rc+ corresponding to 
the maxima in the free energy as a function of T/Tc for temperatures below the critical 
temperature (T < Tc). We see in both the graphs 4c and 4d that the critical radius goes 
down with temperature. In figure 4c we see that the critical bubble radius goes to zero when 
(3 > 1 (long dashed curve having 7 7^ 0). If no curvature term is present then the critical size 
does not go to zero. Even with curvature term, large a does not allow critical radius to go 
to zero (solid curve). Figure 4d plotted for a different set of parameters has similar features 
except that for B^/^ = 300MeV, a ~ IMeV / frn? the critical radius below T^. is zero. The 
only bubbles that can expand are the ones already nucleated above T^.. 

Finally, figures 5a and 5b show the hadron fraction /i as a function of time t as the phase 
transition proceeds to it's completion. Figure 5a shows the overlap for one parameter set 
given there (solid and dashed curves) and clearly shows that if we decrease the surface tension 
the phase transition proceeds slowly (long dashed curve) unlike the case in early Universe 
studies. The effect of 7 has a mixed behavior depending on the value of the bag pressure. 
Figure 5b shows that for low bag pressures, the effect of 7 is to delay the phase transition 
(solid curve) but for larger values of B we see that the curvature correction actually speeds 
up the process of hadronization (long dashed curve). 

The idealized Maxwell construction with the correction to the QGP as well as the the 
hadron gas that we have studied here shows that the time range of the phase transition is 
actually less than what has been given in p|. For the same parameter values as given in 
jHj we find that in the ideal case with corrections, the phase transition ranges from 2.5 fm 
to 22 fm whereas without the corrections it ranges from 3 fm to 37 fm. If we compare our 
study with this corrected ideal Maxwell construction, we find that for S^/^ = 235MeV^ the 
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presence of the curvature term delays the transition to 31 fm for a = bOMeV fm ^, to 53 fm 
for a = 20MeVfm-^ and to more than 60 fm for a ~ 7MeVfm-^ . For B^/^ = 300MeV the 
ideal case streches from 1.9 fm to 6.8 fm but curvature introduces a delay and the transition 
extends to 12 fm for a = bOMeVfm"^, to 16 fm for a = 20MeV fm'"^ and to 27 fm for 
a ~ TMeVfm^'^. It is interesting to note that the ideal Maxwell case range is 1.5 fm to 
2.8 fm for B^/^ = 350MeV, but a ~ 7MeV fm ^ with the curvature term produces a range 
from 1.5 fm to 16 fm. 

4 Conclusions 

In conclusion, we have presented in detail the mechanism of the phase transition of the 
QGP fireball from its creation to it's subsequent hadronization in an ultrarelativistic nuclear 
collision. The inclusion of the interactions in both phases along with the the curvature energy 
term introduces subtlities that are both interesting and quantitative in nature. Broadly, we 
can summarize that whereas interactions lower the critical temperature, the combined effect 
of interaction and curvature energy is to affect both the degree of supercooling and the 
extent of the delay in the phase transition. Smaller values of the surface tension can enhance 
supercooling to almost 50 percent and also increase the delay in the hadronization process 
whereas larger values of the bag pressures speed up the hadronization process, a^^^ = 65MeV 
(a ~ 7MeV"/m~^) is a reasonable value according to the literature and anything above 
a^^^ = lObMeV (a ~ 30MeV fm~'^) seems to be on the higher side |15j. The critical 
temperature is decided by the bag pressure and the inclusion of interactions in the two phases 
suppresses it. Tc ~ 170 — 180Me\^ is favored in the literature, so a bag pressure in excess 
of 300MeV seems to be on the higher side. Inclusion of the curvature effect has a mixed 
behavior in that, for larger values of B the transition is speeded up whereas for smaller values 
of B delay is enhanced. Another interesting observation is that the hadronization actually 
begins the moment the QGP is created, through equilibrium sized bubbles and is later 
joined by the expanding bubbles below the critical temperature. For reasonable parameter 
values the nucleation of bubbles below Tc can get stopped when (3 > 1 and now even the 
equilibrium sized bubbles join the expansion to complete the phase transition. Finally, for 
the ideal Maxwell construction, inclusion of interactions reduces the range of the transition 
but the presence of the curvature term introduces a delay in the completion of the transition 
compared to the ideal case. 
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Figure Captions 

Fig. 1. Free energy W in units of MeV as a function of bubble radius r in fm. In fig. la Sohd 
(with 7) and dashed (without 7) curves are for B^^^ — 235MeV and a — bOMeVfrn'"^. 
Long dashed (with 7) and dotted (without 7) curves are for B^^^ — 235MeV and a — 
20MeVfm-'^. In fig. lb Sohd (with 7) and dashed (without 7) curves are for = 235MeV 
and a ~ 7MeVfm~'^. Long dashed (with 7) and dotted (without 7) curves are for S^/^ = 
SOOMeV and a ~ 7MeVfm~'^. Both the curves la and lb are at a temperature of LOlTp. 
In the lower panel the curves Ic and Id are labelled as in Fig. la and lb respectively except 
that both are at a temperature of .QQTc- 

Fig. 2. Temperature T/Tc as a function of time t in fm. In fig. 2a and 2b Solid, dashed, 
dotted and long dashed curves are labelled as in fig. la and lb respectively. 

Fig. 3. Log of the nucleation time T(Tau) in units of fm/c as a function of temperature. 
The figures 3a and 3b are labelled as in fig. la and lb respectively. 

Fig. 4. Upper panel shows critical bubble radius Vc in fm as a function of T/Tc. They 
are the minimum free energy equilibrium critcal bubbles. The solid and long dashed curves 
are both with 7 correction. These bubbles do not exist without the curvature correction. 
Figures 4a and 4b are labelled as in fig. la and lb respectively. Lower panel shows critical 
bubble radius Tc in fm as a function of T/Tc. They are the maximum free energy expanding 
critcal bubbles. Figures 4c and 4d are labelled as in fig. la and lb respectively. 

Fig. 5. The hadron fraction h as a function of time t in fm. Curves in fig. 5a and 5b are 
labelled as in fig. la and lb respectively. 
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Figure 1: Free energy W in units of MeV as a function of bubble radius r in fm. In 
fig. la Solid (with 7) and dashed (without 7) curves are for B^^"^ — 235MeV and a — 
bOMeVfm"'^. Long dashed (with 7) and dotted (without 7) curves are for B^^^ — 235MeV 
and (7 = 20Mey/m"^. In fig. lb Solid (with 7) and dashed (without 7) curves are for 
51/4 = 235MeV and a ~ TMeVfm-'^. Long dashed (with 7) and dotted (without 7) 
curves are for 5^/^ = SOOMeV and a ~ IMeVfrir^. Both the curves la and lb are at a 
temperature of l.OlTc. In the lower panel the curves Ic and Id are labelled as in Fig. la and 
lb respectively except that both are at a temperature of .99Tc. 
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Figure 2: Temperature T/Tc as a function of time t in fm. In fig. 2a and 2b Solid, dashed, 
dotted and long dashed curves are labelled as in fig. la and lb respectively. 
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Figure 3: Log of the nucleation time T(Tau) in units of fm/c as a function of temperature. 
The figures 3a and 3b are labelled as in fig. la and lb respectively. 
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Figure 4: Upper panel shows critical bubble radius Tc in fm as a function of T /Tc- They 
are the minimum free energy equilibrium critcal bubbles. The solid and long dashed curves 
are both with 7 correction. These bubbles do not exist without the curvature correction. 
Figures 4a and 4b are labelled as in fig. la and lb respectively. Lower panel shows critical 
bubble radius Tc in fm as a function of T /Tc- They are the maximum free energy expanding 
critcal bubbles. Figures 4c and 4d are labelled as in fig. la and lb respectively. 
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Figure 5: The hadron fraction h as a function of time i in fm. Curves in fig. 5 a and 5b are 
labelled as in fig. la and lb respectively. 
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